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Abstract 



We show that with every separable calssical Stackel system of Benenti type on a Riemannian space one 
can associate, by a proper deformation of the metric tensor, a multi-parameter family of non-Hamiltonian 
systems on the same space, sharing the same trajectories and related to the seed system by appropriate 
reciprocal transformations. These system are known as bi-cofactor systems and are integrable in quadra- 
^ ' tures as the seed Hamiltonian system is. We show that with each class of bi-cofactor systems a pair 

fT^ \ of separation curves can be related. We also investigate conditions under which a given flat bi-cofactor 

system can be deformed to a family of geodesically equivalent flat bi-cofactor systems. 

AMS 2000 Subject Classification: 70H06, 70H20, 37J35, 14H70 

^ ■ 1 Introduction 
O 

A significant progress in the geometric separability theory for classical Hamiltonian systems separable by 
Hamilton-Jacobi method has been achieved in recent years (see for example [1] - [4]). Among other things 
a new class of non-Hamiltonian Newton systems was introduced [S], [B]. These systems were shown to have 
$— ( ' very interesting geometric properties when considered as systems on Riemann spaces [I], [5] (see also [5]). In 
[lU] we showed that they can be separated by the Hamilton-Jacobi method after certain reparametrization 
of the evolution parameter (see also jTT]). Originally these systems were called quasi-Lagrangian systems. 
In the present literature they are called bi-cofactor systems or cofactor-pair systems. In [12] it was further 
shown that each bi-cofactor system is geodesically equivalent (in the classical sense of Levi-Civita |13) ) to some 
separable Lagrangian system which means that it has the same trajectories on the underlying configuration 
manifold as the Lagrangian system only traversed with a different speed and moreover that the metric tensors 
associated with both systems are equivalent i.e. have the same geodesies (considered as unparametrized 
curves). In the same paper one can also find a thorough geometric theory of bi-cofactor systems on an 
arbitrary pseudoriemannian space. 

In the present paper we demonstrate on the level of differential equations the geodesic equivalence prop- 
erties of cofactor and bi-cofactor systems expressed by an appropriate class of reciprocal transformations. We 
clarify and systematize their bi-quasihamiltonian formulation on the phase space. We show explicitly that a 
bi-cofactor system is geodesically equivalent to two different separable Hamiltonian systems of Benenti type 
and we show explicitly the transformation between all geometric structures associated with these two Benenti 
systems and the original bi-cofactor system. We further demonstrate that with each bi-cofactor system one 
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can relate two different separation curves and we find a map between tliese curves. From this point of view we 
therefore show that with each pair of separation curves that are related through the above mentioned map we 
can associate a whole class of geodesically equivalent bi-cofactor systems. Every such class contains at least 
two separable Hamiltonian systems and on the phase space all the members of a given class are related by 
a composition of an appropriate non-canonical transformation and a reciprocal transformation. Further, we 
investigate geodesically equivalent families of flat (in the sense of the underlying metric tensor) cofactor sys- 
tems and find a sufficient condition for a so called J-tensor to generate from any given flat bi-cofactor system 
a multi-parameter family of flat bi-cofactor systems. Finally, we illustrate our considerations by presenting 
a thorough example of the class of separable bi-cofactor systems geodesically equivalent to the Henon-Heiles 
system and then specify this example to the flat case. 



2 Cofactor systems 

Let us consider the following Newton system 



= z = l,...,n (1) 

where q^ are some coordinates on an n-dimensional pseudo-Riemannian manifold Q endowed with a metric 
tensor g = {gij) and where F — {F^) is a vector field on Q representing the force which we assume time- 
and velocity-independent. Here and in what follows we use the Einstein summation convention if not stated 
otherwise. The functions F'-j, are the Christoffel symbols of the Levi-Civita connection associated with the 
metric tensor g and if all F*^ are zero we call the system ([T]) a flat Newton system. In case that F = (P) is 
the equation of geodesic motion on Q and the variable t becomes an affine parameter of geodesic lines. 
If the force F is conservative (potential) i.e. if 

F = -VV = -GdV, (2) 

where G = g~^ is the contravariant form of the metric tensor g and where V = V{q) is a potential function, 
then ([1]) is equivalent to the Lagrangian system 

d dC dC i d ^ . 

dtd?-W^'' ' ^dt'^^ ^ = l,...,n (3) 

on the tangent bundle TQ endowed with coordinates {q, v) = (q^, . . . q", v^, . . . , v"), where C — \gij{q)v'^v^ — 
V{q) is a Lagrangian of the system. By the Legendre map pi = gijv'^ the system is transformed to the 
Hamiltonian dynamical system 



d 

dt \ p 




s| \=n,dH (4) 

dp J 



on the cotangent bundle T*Q endowed with coordinates {q,p) — {q'^,Pj) where H — \G^^ {q)piPj + V{q) is the 
Hamiltonian of the system, H^ is the canonical Poisson tensor and dH is the differential of H . 
We will now remind the notion of a J-tensor. 

Definition 1 A (l,l)-tensor J — (Jj) on Q is called a J-tensor associated with the metric g or G = g~^ 
(we often write that 3 is a Jg -tensor when emphasizing the underlying metric) if its contravariant form 
J'-* = JIG''^ is a symmetric {2,0)-tensor and if 3 itself satisfies the following characteristic equation 

VhJ}^{ajSl + a'g,h) (5) 

where V/i is the covariant derivative associated with the metric g and where Ui is some 1-form. 

From ([5]) it follows that the Nijenhuis torsion of J vanishes: 

4V|,i4-j/^V[,4]=.o 

(the square brackets denote skew-symmetric permutations of indices the index h is not permuted) and 
that J is a conformal Killing tensor of trace type which means that Jij = Jigkj satisfies the relation V(ft J^) = 
(^{hdij) with ai = 9itr J (the brackets denote symmetric permutations of indices h,i,j). 
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Remark 2 All J -tensors of a given metric tensor g constitute an R-linear vector space of dimension less or 
equal to ^{n + + 2). This space attains its maximum dimension for metrics of constant curvature. In 
case the metric g is pseudoeuclidean so that g = diag(£i, . . . ,e„) with Si = ±1 in its Cartesian coordinates, 
the general form of J in these coordinates is \12}j 

= mq'q^ + (3'q^ + (3^ q' + 7*^' (6) 

where m^(3'^ and^''^ = 7^' are ^(n+l)(n+2) independent constants and where — JIG^^ is the contravariant 
form of 3 . 

If a J-tensor J has n real and simple eigenvalues then it is called L-tensor and its signed eigenvalues 
(A^, . . . , A") given by det (J + A(q)/) = define a coordinate web on Q. Such webs will turn out to be 
separation webs for our systems (see below). See [T2] for further details on J-tensors and L-tensors. 

The system ([IJ is called cofactor if the force F has the following form 

F = -(cof J)"^ VF (7) 

for some J-tensor J, where cof J is the cofactor matrix of J (i.e. the transposed matrix of signed minors of 
J) so that Jcof J — (cof J) J — (det J) / or in case that J is invertible cof J — (det J) J^^. In the case 3 = 1 
the system ^ becomes Lagrangian (potential). 

In our further considerations the notion of equivalent metric tensors will play an important role. Two metric 
tensors G and G on manifold Q are said to be equivalent if their geodesies locally coincide as unparametrized 
curves. As it was shown in [12 , a metric G admits an equivalent metric G if and only if it admits a nonsingular 

J-tensor J. In such a case G*"' = aJ^G''^ = crJ*-' or in the matrix form 

G = a3G 

with a = det J = ^ where t and t are affine parameters associated with the (parametrized) geodesic of G 
and G respectively. Moreover, J^^ is a J-tensor for the new metric G. 

Two dynamical systems {g, F) and (g, F) of the form ([T]) on Q are said to be equivalent if their trajectories 
coincide up to a reparametrization of the evolution parameter. Moreover, they are called geodesically equivalent 
if also metrics g and g are equivalent. As it was proved in [12] two systems (g, F) and (g, F) are geodesically 
equivalent if and only if the metric g admits a nonsingular J-tensor J such that 

G'^ = aJlG^^ , F = a^F, cr = det J. 

Also in this case the evolution parameters t and t of the systems {g, F) and (5, F) are related through the 
above mentioned reciprocal transformation 

dt 

-= = (T. 

dl 

We will now show that every cofactor system belongs to a whole class of geodesically equivalent cofactor 
systems. 

Theorem 3 Consider the cofactor system 

Assume that Ji is another J-tensor for the metric G and denote by Gi = criJiG (with o\ = det Jij a new 
metric tensor equivalent to G. In a new independent variable ti defined through the reciprocal transformation 

dt 

dti = — 

CTl 

the cofactor system (0^ attains the form 

f -<^'"«'f^-(H("r')l-v"V)-, „ p) 

where (r*-^^)^j, are Christoffel symbols of the metric Gi and V^^^ = Gid is the gradient operator associated 
with the metric Gi. 
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Proof. Since dti = dt/ai we have, by the cham rule, 

dq^ 1 dq^ d'^q^ 1 d'^q^ 1 dg' dai dqi 
dt (Ji dti ' dt"^ dt\ dti dqi dti 

Moreover (see for example [13]) the ChristofFel symbols of G and Gi are related by 
Further 

V^^V = GidV = (TiJiGdV = (JiJiVy, 

so that 

(cof J)"^ VV=— (cof J)"^ JrV^^V = \ (cof J)"^ cof JiV^^V = 
CTi 

= ^ cof(J-i) cof Ji V^^V = ^ cof(JiJ-i)V(iV = 
= i,(cof(JJ-)-^)v«y. 
Plugging all this into dH) we obtain m 

Remark 4 The tensor J is a Jci-tensor i.e. a J -tensor for the metric Gi. It means that the system (0) 
is a cofactor system geodesically equivalent to with Gi as the underlying metric. 

Note that in the particular case Ji — J the system © becomes potential 

d^q' , dqi dq'' .i 

with the affine parameter 



^ dt 

dti = at = — 

(T 

and with F^j, and V — Gd defined by the new metric G = aJG with a = dct(J). This shows that every 
cofactor system is geodesically equivalent (in the sense of the definition above) to a potential system. This 
fact yields us a possibility of determining a quasi-hamiltonian formulation for the cofactor system ([8]). 

Proposition 5 The cofactor system Q) has on T*Q the following quasi- Hamiltonian representation: 



q 

dt \ p I a 



= -TlncdH (12) 



with the noncanonical Poisson operator 

-j^ n J' " [J^f ' 'dq 

and with the Hamiltonian 

H{q,p) ^ \f{coi 3)Gp + V{q). (13) 

Proof. The systems ([5]) and (|lip are related by the reciprocal transformation di. = dt/a with a = a{q) 
yielding that dq^ /di = adq^/dt. Let us thus introduce new variables on TQ: 

q — q, V = av. (14) 
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The Lagrangian of (|lip written in coordinates (5, v) is: 



This Lagrangian defines a new Legendre map from TQ to T*Q that is just the fiberwise isomorphism between 
TQ and T*Q induced by the new metric 'g i.e. p — 'gv. From G = a3G we have 

— -1 1 ,-1 I/tTn 



so that 



Thus, the map p4p on TQ induces the following non-canonical map on T*Q: 

q = q, p= {3^y\. (15) 
In the coordinates {q,p) the system ([TT|) has the following canonical Hamiltonian representation (cf (jl])): 

^ ! ) - n,di7, (16) 



with the usual Hamiltonian H = ^jF Gp + V{Tj). In order to obtain the quasi- hamiltonian form (fT2|) of ([8]) 
it is enough to transform the system (fT6|) back to the variables {q,p,t). The map between these variables is 

q = q, p={J^y%, dt^^- (17) 

or equivalently 

q — q, p — J^p, dt — adi. (18) 

Note that this map consists of a "space" part (fTTj) that involves only {q,p) and {q,p) variables followed by the 
reciprocal transformation (rcparametrization of evolution parameter) dt ~ dt/a. By using that d/dt — ad/dt 
(which generates the factor l/cr in (fT2|) ) and after some calculations that exploit the fact that J is torsionfree 
we obtain ^\2\ with H denoting the function H written in (q,p)-coordinates. Since 

f^Clp = p^J" VJG(J^)" V = P^crJ-^Gp = p^ cof(J)Gp 
we get that H is of the form (fT3|) . ■ 



3 Bi-cofactor systems 

The system of Newton equations of the form 
d'^q^ j dq^ dq^ 



dt^ 



r}.^^ = - ((cof Ji)-^ Vl^) = - ((cof J^)-^ VW^) (19) 



with two independent Jc-tensors Ji and J 2 and with two different potentials V and W is called a bi-cofactor 
system on Q. It means that the force F has two different cofactor representations of the form ([7]). The 
following is a simple corollary of Theorem [3l 

Proposition 6 Assume that the metric G has a third J -tensor J3 and denote by G3 ~ (73J3G (with (T3 = 
det J3J a new metric tensor equivalent to G. In the new independent variable t^ defined through 

dts = - (20) 
CT3 
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the hi- cof actor system ififf|) attains the form 



dt| ^ '^'^ dta dts 



[cot(J2J;')]"' V<»lH')' (21) 

where (r*-'^^)*j. are Christoffel symbols of the metric G3 and V^"^-* — G^d. 

As before, both the tensor Ji Jg"^ and J2 Jj^^ are Jcg-tcnsors so that Gy, is the underlying metric of the 
system (1^ . 

In case that J3 = Ji the system ((2T|) attains the potential-cofactor form 



r,.i5;^ = -(wr.-({cofD-'w)' (22) 



df ' dt dt 

with the affine parameter dt — dt/a and with J = J2 J]^ being a Jg-tcnsor for the new metric G = ctJiG 
with a = det(Ji) = cti. 

If J3 = J2 then the system (PT|) attains the cofactor-potential form 

f-f^'-ff -)'-(™')" 

with the affine parameter dt = dt/a and with J — Ji = J being a Jg-tensor for the new metric 
G = (t32G with cf = det(J2) = CT2. 

Proposition 7 T/ie bi-cofactor system I119\ ) has on T*Q the following bi-quasihamiltonian representation: 

^( l\= ln„,(Ji)di/ = —Unc{J2)dF, (24) 
dt \ p J ai (T2 



with two compatible noncanonical Poisson operators n„c(Ji) and n„c(J2) given by 

and with the Hamiltonians 

H = ip^(cof Ji)Gp + V{q), F = ip^(cof J2)Gp + W{q). 

The representation ([M]) foUows directly from Proposition [S] applied independently to both cofactor rep- 
resentations of (Uni). The fact that the operators n„c(Ji) and I\-nc{^2) are compatible (i.e. that any linear 
combination ?7in„c(Ji) + ri2^nc{^2) is Poisson) is shown below. In the particular case of potential-cofactor 
systems (j22p and (1231) this proposition yields their well-known quasi- bi-Hamiltonian representation [15j , [16j . 

Theorem 8 1. The system \24^ has n constants of motion 



Hr =Er + Vr{q) = \p^ K^Gp + Vr{q), r = 1, n, (25) 

(with H = Hi and F = Hn) where Kr are (1, \)-Killing tensors (for the metric G) defined by 

co/(J2+CJi) = E^»-«^' (26) 

(so that Ki = cof Ji, Kn ~ cof J2j and where the potentials V,- can be obtained from two equivalent 
formulas 

VVr = —KrJlVVi or VVr = —KrJ2'^V„, V ^Vi,W = (27) 

(71 fT2 
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2. The constants are in involution with respect to both operators n„c(Ji) o.'rid n„c(J2)' 

{i7r,i?Jn„,(Ji) {-^■r,i?s}n„,(j2) = Jor all r,s = l,...,n. 
To prove this theorem, we will first need 
Proposition 9 In the variables {q,p,t) related with {q,p,t) through the map 

q^q, p^{3^)-\, dt^^ (28) 

the system {24\j attains the quasi-bihamiltonian form 

with H = H and F = F (as functions on T*Q) and with 

n, = n„e(Ji), n„,(j) = n„e(J2) (30) 

(as tensors on T*Q). Moreover, the tensor He is canonical in {Ti,p) -variables. Similarly, in the variables 
(g, p, t) defined by 

q = q, p= {3^y\, dt = —. (31) 

attains the form 

1 



-Unci3)dH ^n^dF, (32) 



dt \ P J dct(J) 

with H = H and F = F (as functions on T*Q) and with 

n, = n„,(j2), n„c(J) = n„,(Ji) (33) 

(again considered as tensors on T*Q) so that 

nc = n„,(j), ne = n„e(J). 

Again, the tensor lie is canonical in {q,p) -variables. 

This proposition can be proved either by direct calculation or by observing that the underlying bi-cofactor 
system (fT9|) has in the variables {q, t) the potential-cofactor form ([22|) and in the variables (q, t) the cofactor- 
potential form and using arguments similar to those used in the proof of Proposition [S] 
Proof, (of Theorem [8]). By Proposition [9l the system p4l) has in variables {q,p,t) the form ((29)) so that it is 
a so called Benenti system and therefore (see [3]) He = n„c(Ji), n„c(J) — n„c(J2) are compatible and the 
system has n constants of motion of the form 

Hr = Er + Vriq) = \p^KrGp + F,(g), r = 1, n, (34) 

with G = (TiJiG and where the Killing tensors Kr of the metric G are determined by the expansion 

n-l 

cof(j+e/) -^i?„-.r (35) 

(so that Ki = I, Kn — cof J) while Vr are separable potentials satisfying 

Kr^Vi^^Vr- (36) 

On the other hand, 
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cof (j+e/) = cof (JaJr'+e/) = cof(Ji)-i cof (j2+eJi) = i^r' 

SO that by comparing with (|35p we obtain 

K, = K-^K„ z = (37) 

and 

The last equahty follows from J^^K^^ — ail and JiG (J^^)"^ = JiJ^^G = G, so that indeed Er = fi'r if 
we define Ki as in ([^^ . Thus, if we put V^r = V^^ we obtain that iJ,. = i/^ (as functions on T*Q). Now, 
substituting (|57)) into we get 

i^f^if^CTi JiVV^i = (Ti JiVK or KrJiVVi = XiJiVK 



which yields the first formula in ([27| . Naturally, the functions iJ^ Poisson-commute with respect to both Pois- 
son tensors n„c(Ji) in (|24| since i/^ = -ffr Poisson-commute with respect to both Poisson tensors lie, n„c(J) 
in ((29|) and since these tensors are just n„c(Ji) written in the variables {q,p), according to ([30|) . That proves 
all the statements in Theorem [8] except the second formula in ([27l) . Consider now the system ([32]) . It is also 
a Benenti system so it also has n constants of motion of the form 

Hr=Er+ Vr{q) = ^p^ KrGp + Vriq), T = 1, n, (38) 

(with G = (T2J2G and with H = F = Hi) where the Killing tensors of the metric G are determined 

by 



cof (j+c/) = J2 



(so that Ki = J, = cof J) while Vr are separable potentials satisfying 

KrVVi = yVr. (39) 
Repeating the procedure above we obtain an equivalent proof of Theorem [51 However, this time we get 

n-l 

Kn-^e = cof (j+e/) = cof (Ji + = cof (Ji + ej2) = 

i=0 

n—1 n— 1 

which gives ^ 

K, = K-^K^^^+i, z = l,...,7i. (40) 

This also yields Ei = En-i+i and thus Vi = V^-j+i, = Hn-i+i = H„-i+i for all i = l,...,n. By 
transforming the formula (j39[) to ((7,p)-coordinates (similarly as we did for bar-coordinates) we obtain the 
second formula in ([27)l . Finally, the functions Ki defined by ([26)) must be (1, 1)-Killing tensors for G since 
cof(J2+^Ji) is a (1, 1)-Killing tensor for any value of the parameter ^ and since Killing tensors of G constitute 
a vector space. ■ 

Note also that the direct map between variables (g, p, t) and (q, p, t) is obtained by composing the map 
with the map ([?T|) . It attains the form 



g, p= J p, rft = or g = (7, p= J p, d? ^ ~ ■ (41) 

^ ^ det(J) V / dct(J) 
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Further, by comparing p7|) and pOj) we obtain 

In the remaining part of this chapter wiU shortly discuss how the two equivalent systems (|29p and (|32p can 
be embedded in quasi-bihaniiltonian chains and discuss the relation between these chains. 

Since the system (j29|) has n commuting with respect to both operators He and n„c integrals of motion 
it belongs to the set of n commuting Hamiltonian vector fields 



— d 



UcdHr = Xr, r^l,...,n, (42) 



(where dii — dt — d</det(Ji)) and the system (|29p itself defines the first vector field Xi. Similarly, since 
the system ([5^ has n commuting with respect to both He and n„c integrals of motion Hr it belongs to the 
set of n commuting Hamiltonian vector fields 



dtr V P 



:=:UcdHr=Xr, r = l,...,n, (43) 



(where dti = dt = dt/det(J2) so that dti — c?<i/det(J)) and it also is the first vector field Xi. By the above 
construction, the vector fields Xi and Xi are parallel 

Xi = det(J) Xi or Xi = det(J) Xi 

which once again refiects the geodesic equivalence of the systems (f22|) and (f23|) on Q. 
Moreover, vector fields (|42p belong to the following quasi-bihaniiltonian chain: 

Xi = n, dHi = ^n„,(j)di?„ 

Pn 

Xr - n, dHr - '^Tlnc(S)dHn - Tlnc{J}dHr-l, T = 2, .., n, 
Pn 

where the functions are defined through the polynomial expansion of det(J + ^/): 

n 

det(J + a) = ^p,,r"^ 



(so that p„ = det J). Similarly, vector fields (|43p belong to a similar quasi-bihamiltonian chain: 

= u,dHi ^ ^n„,(j)di?„ 

Xr = UcdHr = ^n„,(J)di/„ - Unc(3)dHr-l, r = 2, ...,n, 
Pn 

where pr are defined through 

n 

det(J + e/)=^P.r-^ 

4=0 

(so that Pn — det J). Since J = J ^ we have that and pr are related via 

~ Pn — r — Pn — r 

Pr = — — or p,, = — 

Pn Pn 

Comparing both chains we obtain that the vector fields Xr and Xr are related through 

^1 = Pn ^» — Pn-i+l^l ^ -'^n-i+2, i = 2, . . . , 71. 
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4 Flat bi-cofactor systems. 



Let us recall that a pseudoriemannian space is called space of constant curvature if the curvature tensor Rijki 
has the form 

Rijki = K {gjigik - gjkgu) (44) 

for some scalar function K . By Bianchi identity it follows then that if is a constant, related to scalar (Ricci) 
curvature x = Ritg^^ through k = Kn{n — 1). Thus, for such spaces the condition k = ox K = implies 
that the Riemann tensor Rijki is zero i.e. that the metric g is flat. 

Suppose now that g is a metric of constant curvature and that 'g is another metric tensor obtained by 
deforming g through 

G = aJG (45) 

(with J being a Jc-tensor J and with a = det J). Then, by the classical result of Beltrami |17| we know 
that g is also of constant curvature. Moreover, for two metrics g and g that are geodesically equivalent and 
of constant curvature their scalar curvatures >c and 35 are related by the formula 

= X9ij - ^ifi + fifj (46) 
(see [H] p. 293) where the covector fi is defined as 



2(n+ 1) % V detg 
A simple calculation shows that for our choice of g,^ we have 

Ji = —Oi where = - — . 

n + 1 fj oqi 

Substituting this into (|46p and performing contraction with G we obtain 

><tr J + -^J""^ (o-iCTi + -A-rV^a. 

n+l \ ^ J ' n+1 * J 



_ a 

>c — — 

n 



(the summation convention applies as usual). Thus, we see that if >f = then a sufficient condition for x to 
be zero IS 

J'' (^^^. + l^i^^'^j) = 0- (48) 

Let us now assume that the metric G of the system is flat (i.e. >f = 0) so that in some coordinate system 
((?') it assumes the form 

G = diag(ei, . . . , e„) with ~ ±1. (49) 

(note that then g = G^^ — diag(ei, . . . , £„) in this particular coordinate system too while F*^ = 0). Suppose 
now that we want to "deform" this system as in Proposition [6] by introducing the new independent variable 
dts = dt/a where a = det J for some new J^-tensor J but in such a way that the resulting equivalent metric 
G = aJG is also flat so that the geodesically equivalent system (|2T|) is a flat Newton system (in this section 
we will use J, a and G instead of J3, det J3 and G3 to shorten the notation). A sufficient condition for doing 
this is to take J that satisfies In the Cartesian (with respect to g) coordinates (q*) the contravariant 

form of the tensor J is given by ([6]). However, by Theorem B.4.3 in [12] we know that m = or else 35 7^ 0. 
Thus, our aim is to find a more explicit form of the condition (I48|) for J given by Q with m = 0. 

Let us for the moment denote the (2, 0)-form of J as given in ^ by Jc (J-contravariant) so that Jc = JG 
or J = Jcg- We have then 

Theorem 10 Assume that G is of the form J^ffp and that J is a jQ-tensor such that its contravariant form 
Jc is given by (0). Then for the metric G — aJG to be flat it is sufficient that m = and 

^^g(cof J) /3 = or /3^(cofJc)/3 = 0. (50) 
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Proof. Both conditions in (|50p are equivalent since cof Jc = cof(JG) = cof G cof J = dct{G) g cof J. We 
have to show that the condition (|48|l in our setting attains the form ([50|) . Since J is torsionless it satisfies the 
identity 

9(trJ) ,,da 



or in the matrix form 



crdftr J) ^ J^dcr. 



(51) 



Since Jj = J guj = P^£j<l-' + P-'^jQ^ + 7'"'ej (no summation) we have tr(J) = — 2P^Siq^ + Y^Si so that 
d(tr J) = {si(3\ . . . ,£„/3")^. Thus, jH]) reads 



Therefore 



Further 



da = 2{cof3y g(3 = 2<?(cof J)^. 



1 da j,^^da^ ^ X{daf 3 da = ^P^g (cof J) /3. 



■' fj^ dq^ " dqi fJ^ 
Oq^ \ a oqi 



(52) 
(53) 



-aiJ^Ja, + -J'J- 



a dq^dq^ 



But, using ([5^ twice and we obtain 



|-(4(cofj);)-(cofj);|-4 

4^"^^ (cof J) /3 - 2(n + l)/3^5 (cof J) /3. 



13' = 2-^/3^ - 2(n + 1)13'' Ek (cof J)^ 



so that 



Thus, 



J^3 



d'^a 
dq^dq^ 



= 2(l-n)^^g(cof J)/3. 



r^Wia, = 

^ a 



(1 +n)/3^g(cof J)/3. 



(54) 



Plugging ([531) and ([54]) into (gll) we immediately obtain ([501). ■ 

Therefore, we have showed that for any flat bi-cofactor system ([T9l) there exists a multi-parameter family 
(with in(n + 3) — 1 parameters) of geodesically equivalent (but algebraically very different) flat bi-cofactor 
systems. 



(55) 



Remark 11 The condition i50\) can be written as 

f3^g (cof 75) /3 = or (cof 7) f3 = 0. 



5 Separation curves for bi-cofactor systems 

A system of n algebraic equations of the form 

((5j(A*, ^i; ai, a„) = 0, i = l,...,n, det 



doi 



^0, 



(56) 



each containing only one pair (A*,/ii) of coordinates (A,/i) on T*Q (and with real coefficients Oi) is called 
separation relations. The condition in (|56p means that we can solve the equations (j56p with respect to 
obtaining n independent functions on T*Q of the form = Hi{\, fi), j = 1, . . . , n. If the functions Wi(A*, a) 
are solutions of a system of n decoupled ODE's 



dW,{X\a) 
<Pj: 1 A = ,ai, , 



a„ = 0, i = 1, ...,n. 



(57) 
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then the function W{X,a) — X]"=i^j('^N ^ solution of aU the equations ([57)1 and simultaneously it is an 
additively separable solution of all Hamilton- Jacobi equations 



H,{\ ,...,A 1 _ 



i^l,...,n. (58) 



simply because solving ()56|) to the form a.i ~ Hi{X,ii) is a purely algebraic operation. The Hamiltonians Hi 
Poisson-commute by the classical theorem of Jacobi. The function W{X, a) is a generating function for the 
canonical transformation (A, /i) — s- (6, a) to the new set of coordinates that simultaneously linearize all the 
Hamiltonian equations 

ut, =IlcdH., = Xh,, i = l,...,n. (59) 

The coordinates (A,/i) are thus called the separation coordinates for the Hamiltonian systems (j59p . 

In the case that the relations are affine in the obtained systems belong to the well-known class of 
(generalized) Stackel separable systems. 

Let us now consider a special subclass of Stackel systems given by the following separation relations: 

H,{XT-' - H2{XT~' + - + i-ir-'H^ = lf^iy)^^l + 7.(A0, « - 1, -.n, (60) 

where fi and ji are smooth functions. Such systems are known as Benenti systems. In the particular case that 
/i(A') = /(A') and 7i(A') = 7(A*), separation relations ([50)1 are given by n copies of the so called separation 
curve 

i7lA"-l - i/2A"-2 + ... + (-l)"-li/„ ^ ^/(A)^^ ^ ^^^y (gl) 

so that now A and /i G R. By solving the system of n copies of this relation (with i-th copy containing 
variables labelled {Xi,fii)) with respect to Hi we find that the Hamiltonians Hi attain the form 

Hr =Er + Vr{X) = ^^l'^KrG^l + ^(A), r 1, ...,n, 

(cf (|25p) with the metric tensor 

^ ffi^') /(A") 

G = d:ag(^^,...,^ 

where = (A* — A-') while the (1, l)-tensors Kr are generated by the expansion 

ra-l 



cof (J+^/) - J2 



dpr dpr 



1=0 

with the Jc-tensor J — diag(Ai, . . . , A„). Thus [12]: 

Kr = diai 

where the functions pr can be obtained from 

n 

det(J + e/)-^P.r"^ 

i=0 

Coordinate-free expression for Kr is as follows (TB] 

Kr+i = PrI - 3Kr, r = 0, 1 , . . . , n - 1 , po = 1, Ko = 0, 

or alternatively 

r-l 

Kr = ^Pk i-jy^^^'' , r = l,...,n. 



k=0 
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It is important to stress that in case that eigenvalues of J are not simple the obtained tensors Kj. will not be 
independent and thus will not generate an integrable system (see also below). 

For a particular choice 7(A') = (A*)*^, /c S Z in the separation curve ([6T|) we obtain a family of separable 
potentials that can be constructed recursively by [19] 



= pX^^ - x/W with = (-i)"-M, 



This recursion can be reversed 



In both cases we put Vi- = for r < or r > n. The above recursion can be written in a matrix form as 



T^(fc-l) ^ /^^(fc) „ y{k) 
Pn 



(62) 
(63) 



v'^''\\) = r''{\)V^°\ kez, 



(64) 



where t/('=)(A) = (^^f ^ (A), Vi''^ (A))^, V^"^ = (0, 0, (-1)"-^)^ and 



4 k) 



R 



( piW -1 

P2(A) -1 

Pn-l(A) 

V Pn(A) 



\ 





-1 
/ 



(65) 



This recursion is equivalent to (|62|) and ()63|) and is invariant with respect to any point change of variables on 
Q as i? in (j65p is expressed by coefficients of the characteristic polynomial of J. The first nontrivial potentials 



in the positive hierarchy are ^^"''(A) = Pr(A), while for the negative hierarchy Vr 



(-1) 



Ar-i(A)/p„(A). 



Let us now once again consider our systems (|22p and (|23|) and their Hamiltonian formulations (|42|) and (|43|) 
respectively. From now on we will additionally assume that the tensor J (and hence J) has all its eigenvalues 
real and simple (i.e. both are so called L-tensors 12 ). Then the tensors Kr (and Kr likewise) are independent 
and thus (|42|) and (|43)) are integrable. Moreover, both systems belong to the class of separable (in the sense 
of Hamilton- Jacobi theory) systems called Benenti systems. It is known that all Hamiltonian flows in l|42p are 
separable in variables (A,/7) where the new coordinates A* are obtained from the characteristic equation of J: 



det(J + A/) = 



(66) 



(i.e. are (signed) eigenvalues of J) while the corresponding momenta are obtained from — ($' ^)^p 
where is the Jacobi matrix of the map $ : q A given by (l66|) . Similarly, all the flows in pS]) are separable 
in variables (A, /I) where A* are obtained from 



det(J + A/) = 



(67) 



with the corresponding momenta /i* obtained hy ^ — (5*' ^)^p where Vf' is the Jacobi matrix of the map 
: q^X given by ((67l) . 



Theorem 12 The separation variables (A, /i) of 114 and the separation variables (A,/i) of 114 3\ l are related 
by the transformation 

~ 1 _ -i_ 

A* = fii — —A (no summation), i — 1, . . . ,n. (68) 

A 

Proof. By comparing (|66p with (|67p we obtain that A = j-- The map between momenta Jl and JI can be 
found in the following way. We know (cf. ([28]) and ((3T|) ) that p — (Jj) ^ P and p = (J|") ^ p which yields 
p= (j^) ^p. Thus 

li^i^f'-^p^ (*'-i)^(J^)-i$'^7Z^$'J"\*'-i)V 
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Since A = 1/A' we see that ^' = 6$' where 

Q = — diae 



,(A )2 (A )V 

so that ^i'^^ = —$'-1 diag ^(A^)^, . . . , (A")'^^ . Inserting it in the above formula yields 



/I=-(cf>'j-^$'-idiag((A')2,...,(A"f) 



T 



But $'J $'^1 = diag(l/A ,...,1/A ) since it is the inverse of the L-tensor J written in its separation 
coordinates A. Inserting it into the above formula we get the map between momenta as in (|68p . ■ 

According to the remarks above, the Benenti system (1^^ in variables (A,/!) has the separation curve 

Hif'-H2X"'' + ... + i-ir-'Hn = ^7(X)1I' +j(X). (69) 
Similarly, the separation curve for the Benenti system p2p is 

iIiA"-i - H^X"-^ + ... + {-ir-'H,, = ^f(x)J,^+j(x). (70) 



Applying the map ([68]) to the separation curve l|70|) . using that Hr — Hn~r+i and comparing the result with 
((551) obtain 



+ (-l)"-i7(A ')A" \ 



Corollary 13 // the functions /, / and 7,7 satisfy the conditions: 

7(e)-(-i)"-7(r')r+\ m = i-ir^'iiOC, e^R (71) 

then the separation curves I169\) and |70[ ) generate two geodesically equivalent systems of Benenti type parametrized 

by two different evolution parameters t and t such that dt — dt/ det{a), where 'a — pn = nr=iA • The corre- 
sponding families of separable potentials for both systems are related by 

y^''^ = (-1)"-^ or y^'^) = for all r e Z. 

It is known that the metric 

^ nlyif7(A) = ES 
ately that the equivalent metric 



is of constant curvature if and only if /(A) = X]fc=o '^'^'A for some constants Ck- From (|7ip it follows immedi- 



G = diag 

is also of constant curvature, as in this case 



/(Ai) /(A") 
. Ai A„ . 



n+l n+1 



/(A) = (-l)"-V(A-i)A"+i = (-1)"-! ^ cfcA"-^-+i = (-1)"-! ^ c„_fc+iA* 



fe=0 fc=0 
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6 Example: flat bi-cofactor systems geodesicafly equivalent to Henon- 
Heiles system 

Let us illustrate the ideas of this paper on the example of the integrable case of the Henon-Heiles system. It 
has the potential-cofactor form (f22|) : 



^ = = - (cof J)-^ W = -^^^<^) _t_| ) j (72) 

with G = / (so that r*j, — and coordinates q are Euclidean) and with the Jg-tcnsor J of the form 

J 



-W 



The potentials V and W are 



The system ([7^ has the quasi-bihamiltonian representation ([M)) 

with Hamiltonians Hr of the form (p4|) . Explicitly: 

Hi=H=\ {Pif + \ fe)' + V{q), H2 = F= \tp{P2 - \q\p2? + W{li). 
The non-canonical Poisson operator n„c reads explicitly as 

The system (|73p separates in variables (A, Ji) that can be found from the characteristic equation (|66p and are 
given by 



7-1 ^2 



Pi = - ^1 =2— h P2 = V ^ ( =1 =2 

yA — A/ VA — A 

while the separation curve (|69p generating Hamiltonians i/^ is 

i?iA - i?2 = ^AJI^ - A^. (74) 

Let us now take another, arbitrary J^-tensor J3. Since the metric G — I in variables ijf' ^'(f') the most 
general form of J3 is ^ which reads now as 

' 1^ mtt + Pit + P2t + 1'' m{t)' + 2/329' + 1^' J 

with 6 arbitrary constants m,Pi,j''^ = 7-'*. Using Proposition [5] (with Ji = / and J2 — J) we see that in a 
new independent variable defined by 

dt 



det (J3) 
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our potential-cofactor system (jTS]) attains the bi-cofactor form pT|) with (r'^'^))^^, being Christoffel symbols of 
the metric G3 — (det J3) J3G. They can be obtained from (fTO|) which reads now (since F^j. = 0) 

rr(3)v' 1 / , 9(dctJ3) ,9(detJ3) 



2dctJ3 V"^ 99fc ' d^j 

It is important to stress that for aU the choices of J3 the obtained system has on Q exactly the same trajectories 
as Henon-Heiles system, only traversed with different speed. Moreover, the metric G3 is of constant curvature 
since it is geodesically equivalent to the flat metric G = I [IT]. 

Among all possible choices of the deforming tensor J3 there is only one that leads to a new potential- 
cofactor system, namely J3 = J. This choice leads to cofactor-potential system with J = J and with 
the metric G — (det j) J G. The metric G is flat since the deforming tensor J3 = J satisfies the conditions 
(jSU]) and Explicitly, we have 

(where of course = g*). This system has the quasi-bihamiltonian form (|32|) with the Hamiltonians as in 
(j38[) {Hi — H2 and H2 — Hi ) where the new momenta p are related with the old momenta through the map 
([4T|) and read explicitly as 

Pi = -q^pi - 2^?-^' P2 = -~2^Pi- 

Our new system separates in variables (A, /I) that can be found from the characteristic equation (|67[) and are 
given by 



1 = -[-- + ~7 , q = 



Ai A2y V^AU^ 



p. = -M, (a^Mi - A^M.) , P. - ((Af Mr - (Af /I.) 



Ai - A2 

Our system can be obtained from the separation curve of the form (|70p that explicitly reads as 

i/iA-i?2 = -^A2^2^A-3 

and can also be obtained from the separation curve (17^ by the transformation ([711) . Let us now introduce a 
new coordinates (r^,r^) on Q defined through 

= -2^' 5^ = (76) 
(see [20J). In (r^,r^) the metric G attains the antidiagonal form 

G(r)=(; J) (77) 

(so that (r^,r^) are flat coordinates for G and (r(r))*j, = 0) while the J^-tensor J becomes 

Our cofactor-potential system (geodesically equivalent to ([7^ ) attains in variables {r^^r"^) the flat Newton 
form 
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with potentials 

Viir)=Vir) = 16 ^ \ V^ir) = W {r) = -^^ ^ , 

while the corresponding Hamiltonians are 

Hi{r, s)^H2{r,s) = S1S2 + Vi(r), 

where the momenta (si, S2) are obtained from the point transformation (|76p and are 

2 ^ ^ 4 

Let us finally make another choice of the deforming tensor J3, namely m = 0, 

/31 = 7" = ^ so that 

5gi 2bt + a)' ^3 - (det J3 j J3G - _^3(^i)3 _52(gi)2(26g2 + a) 

The metric G3 is again flat since the deforming tensor J3 satisfies the conditions ([50)1 and ([55]) . The respective 
-tensors for the related bi-cofactor system are 

Let us now perform a parameter-dependent change of variables to the coordinates {x^,x'^) on Q defined 
through 

^ _ 1 2bq'^ + a 2 _ 1 



2 52gi ' 62^1- 

In the metric G3 attains the antidiagonal form ([77]) . so that (a;^,a;^) are flat coordinates for G3 and 

the J(33 -tensors J become 

Hence, our two-parameter family of flat bi-cofactor systems attains in variables {x^ , x'^) the flat Newton form 
(f ( x^ 



dtl \ x" 



(cof Ji)"' v^-'^vy = - ((cof J2)"' v(3)w^y 



with the potentials 



V 



I Xx^^iyX^ f - ^ax^x'^ +a^(x? f\ 1 4x^ + ax^ \ 

8 54(2.2)5 + 2 66(a;2)5 

14(a;i)2-4axia;2 + a2(a;2)2 3 



V 



8 64(x2)4 

14(a;i)2 ~4aa;ia;2 + a2(^2)2 



66(a;2)4 



(79) 



64(a;2)3 



66(^2) 



2\3 ' 



_ 1 (aa;2 - 2a;^)[4(xi)2 - 4aa;ia;2 -t-a2(a;2)2] 1 aa;2 - 2a;i 

- 64(a;2)4 + 16 66(2.2)4 ' 

All the flat bi-cofactor systems in the 2-parameter family ([75)1 are geodesically equivalent to both the Henon- 
Heiles system (|7^ and the Hamiltonian system (|78p. They also belong to the whole family of such systems 
generated by (|75|) . 
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